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Abstract 



A non-linear diffusion equation is derived by taking into account hopping 
rates depending on the occupation of next neighbouring sites. There appears 
additonal repulsive and attractive forces leading to a changed local mobiltiy. 
The stationary and the time dependent behaviour of the system are studied 
based upon the master equation approach. Different to conventional diffusion 
it results a time dependent bump the position of which increases with time 
described by an anomalous diffusion exponent. The fractal dimension of this 
random walk is exclusively determined by the space dimension. The applica- 
bilty of the model to descibe glasses is discussed. 
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There has been a lot of effort in understanding strongly interacting liquids, especially glass- 
forming, undercooled liquids3S. Obviously, the liquid-glass transition is primarily dynamic 
in origin characterized inevitably by a high cooperativity of local processes!. Due to the 
seminal papers of Gotze et ali and independently of Leutheuseri the mode coupling theory 
for dense liquids offers a theoretical tool to study the slow dynamics in the vicinity of the 
glass transition. The mode coupling theory implies that the (ergodic) liquid state can be 
characterized by the decay of correlations whereas below a critical temperature T c predicted 
by the approach it results a crossover to a solid-like behaviour where fluctuations are effec- 
tively frozen in. Furthermore, if additional thermal activated diffusion processes are taken 
into account the correlations decay but on a longer time scaleH§. The extended mode cou- 
pling theoryi'0 leads to a more realistic behaviour of the correlation function. 
The present paper is motivated by another kind of approach to supercooled liquids and 
glasses. It is based on a proposal due to Fredrickson-Andersen (FA)§ to map the sys- 
tem onto an Ising model where the spin variables are related to local regions with a lower 
or a higher density than an averaged one. Such a facilitated model is analyzed above all 
numerically0 but recently also analytically^'^. Different to the conventional kinetic Ising 
model the dynamical rules for spin-flip processes are changed in order to take into account 
the enhanced immobility of the particles due to the gradual freezing with decreasing tem- 
perature. These restrictions can be easily formulated in a Fock space representation of the 
master equation by assuming flip rates which are dependent on the occupation of neigh- 
bouring sites promoting or preventing the processes under consideration. In this manner 
the inherent cooperativity manifested in a cooperative rearrangement of certain regions in 
order to change mobile states within the subsequent time steps into immobile ones and vice 
versa. The inclusion of mobile vacancies in the facilitated modeffi yields to the occurence 
of nearly nonactivated relaxation additional to the main a-process. 

Here, we study another kind of secondary processes. They are due to the ability of the par- 
ticles to diffuse when the system is above the so called glass temperature T g . The diffusion 
constant goes to zero when the temperature is further decreased. Obviously, the mobility of 
the system is more and more reduced whenever the temperature reaches the glass tempera- 
ture. Here, we consider diffusion for T > T g but take into account restrictions which lead to 
the above mentioned diminishing mobility. To this aim particles can jump from lattice site 
to lattice site which gives rise to diffusion on a large scale. To simulate the freezing process 
the hopping rate depends on the occupation of the nearest neighbour site of a given pair of 
sites between which a particle can jump. In such a manner the diffusion is more and more 
confined through the environment. 

To describe this situation we use the mentioned Fock-space representation for the master 
equation under the influence of an exclusive dynamics. Such processes had been discussed in 
different contexts, compare for instance0~il. As already mentioned the starting point is the 
master equation which describes the stochastic process. A systematic approach is achieved 
by a mapping of the stochastic equation onto an evolution operator entering as the infinites- 
imal time evolution operator. This operator can be expressed in terms of second quantized 
operators. Physically it means the considered system is divided in small cells of size 1 (lattice 
points). Every cell is occupied by particles underlying dynamical rules depending on the 
situation in mind. From a quantum point of view two limiting cases can be realized, an 
unrestricted occupation number and a restricted one obeying the Pauli-principle with the 
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occupation numbers and 1. This situation corresponds to a lattice gas representation with 
empty and single occupied sites which can be also interpreted as the two orientations within 
a spin model. 

The case of unrestricted occupation had been discussed firstly by Doill and had been fur- 
ther elaborated inl^ES. Here, we consider the case of a restricted occupation cell. Hence, 
the problem is to formulate the dynamics in such a way that this restriction is taken into 
account. The method offers a wide range of application, for a recent summary see0. 

The starting point is the Master equation which can be derived on quite general grounds. 
It can be written in the form 

d t P(n,t)=L'P(n,t) (1) 

where P is the probability that a certain configuration characterized by the vector ft at 
time t is realized. The evolution operator L' will be specified by the dynamics of the model. 
Followingl!l'0i!~@ the probability distribution P(ft, t) can be related to a state vector | F(t)) 
in a Fock-space according to P(ft, t) = {ft \ F(t)) with the basisvectors | ft). As a consequence 
the Master equation ([]]) can be transformed to an equivalent equation in a Fock-space 

d t | F(f)> = L | F{t)) (2) 

The operator V in (|1|) is mapped onto the operator L given in a second quantized form. 
Up to now the procedure is independent on the used operators. Usually L is expressed in 
terms of creation and annhiliation operators satisfying Bose commutation rulesc2rH3. To 
avoid double occupancy as for instance in a forest fire model!! or a model simulating traffic 
jamil the method can be extended to the case of restricted occupation numbers per lattice 
sitesJlHO introducing Pauli-operators. These operators commute for different points and 
anticommute at the same lattice point. A further extension to a p-fold occupation has been 
discussed recently"^ 

The relation between the quantum formalism and the probability approach based upon 
the Master equation can be found by expanding the vector | F(t)) with respect to the 
basisvectors of the Fock-space 

\F(t)) = J2P(ft,t)\ft) (3) 

rii 

As it was shown firstly by Doiil the average of an arbitrary physical quantity B{ft) is given 
by the average of the corresponding operator B(t) 

(B(t)) = J2P(n,t)B(n) = {s \ B \ F(t)) (4) 

rii 

This rule remains also valid as well as in case of Pauli-operatorsEl and parafermi-operators 
however with different meaning of the state vector | s) introduced in the last equation. 
Remark that the normalization condition for the state function is manifested in the relation 
{s | F(t)) = 1 and the sum-vector {s \ can be expressed by (s |= 5_)(n |. The evolution 
equation for an operator B can be written 

d t {B) = {s\[B,L)\F(t)) (5) 
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To derive the last relation we have used (s \ L — 0. Note that dynamical equations of the 
classical problem are determined by the commutation rules of the underlying operators and 
the form of the evolution operator. 

Now let us demonstrate the applicability of the method to describe diffusive motion under 
the influence of a constraint. In the second quantized formulation the hopping process is 
simply descibed by the evolution operator 

L = Yl J iJ ( d \ d i ~ ( X ~ n i) n i) ( 6 ) 

<ij> 

where the hopping rate J^- between the adjacent sites i and j can also depend on the 
occupation number of the next nearest sites, see below. Furthermore we have introduced 
the occupation number operator 

Because the original problem is a classical one the operators di and d\ satify the commutation 
rules for Pauli-operators. They commute at different lattice sites and anticommute at the 
same sites. 

[di,d]) = 8ij(l - 2n;) 

Using eq.@ and the commutation relations we get the evolution equation for the averaged 
occupation number 

d t (n r ) = ^2(J rj (nj-n r )) (7) 

r(j) 

The evolution operator allows diffusive processes between nearest neighbour sites when- 
ever a double occupancy is avoided. Now let us generalize the model by including the above 
mentioned topological restrictions due to the freezing of the system. To this aim we have 
to include that the hopping rates J r j depend on the occupation number n. As the simplest 
nontrivial assumption we chose 

J rj = J Xrjm (8) 
i 

Here, Xrji is only nonzero when the triple of lattice indices rjl denotes nearest neighbours. 
Introducing the abbreviation 0y = 1 when i and j are nearest neighbours and zero otherwise 
one finds 

Xrjl = 6rj@w(l - Sji) 

Let us visualize the proceeses for the one dimensional case. A particle can jump from site 
i + 1 to the empty site i when the site « — 1 is occupied. On the other hand a particle jumps 
from site i to the empty site % + 1 under the condition that the site i — 1 is occupied. In such 
a manner the hopping process is performed under the influence of attractive and repulsive 
forces as described before. 

The kinetic equation for the averaged density reads 

d t (n r ) = Xrjlinjni - n r ni) (9) 
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There occurs a hierarchy of evolution equations which cannot be solved in general. Here 
the higher order terms are decoupled due to a mean field approximation. In the resulting 
equation we perform the continuum limit for the field n(x,t) = (rij(t)). The corresponding 
evolution equation is expanded upto the order of I 2 , where I is the cell size introduced before. 
The corresponding evolution equation in this order and in a low density expansion is written 
in the form 

J- l d t n = (z- l)nV 2 n - 2(Vn) 2 (10) 

where z is the coordination number and for simplicity we have set / = 1. In case of a simple 
cubic lattice z is simply related to the space dimension via z = 2d. Because the hopping 
rate depends on the concentration itself the dynamical equation is not conserved. 
The stationary solution n s (x) can be found from eq.fllOD 

| Vn s (x) \=n s (x) 2/{z - l) (11) 

For simplicity let us study a radial symmetric solution initiated by a radial symmetric initial 
condition. The exact result is 
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where c and c\ are constants. The special case d = 2 is incorporated since (2 — d) V 2 d 
leads to a logarithmic term 

n s( r ) = (ci + clnr) 3 (13) 

For large distances the solution behaves like 

a -i (z-4)(z-l) (d — 2)(2d— 1) 

nJr) ~ cr~ a with a = - . - . 1 = - ^ (14) 

K J 2(z-3) 2d -3 K J 

where for d — 2 it results in n s oc In 3 r which is different from conventional diffusion n cx lnr. 
Remark that the stationary solution in case of usual diffusion is described by an exponent 
a = d — 2. Therefore d = 2 is a special case which will be also manifested in the dynamical 
analysis. 

Now let us discuss the stability of the stationary solution against a time dependent per- 
turbation. To this aim we write n(x,t) = n s (x) + n\(x,t). (From now on we include the 
coupling constant J introduced in eq. (§) in the time variable t.) It is easy to check that 
the correction term n\ obeys 

d t m = (z-l) (niV 2 n s + n s VV) - 4Vn s Vni (15) 

Such a partial differential equation can be transformed into a diffusion like equation making 
the ansatz n\{x,t) = (f>{x)ip{x,t) where the function <f>{x) is chosen in such a manner that 
the linear gradient term disappears: 

(z - l)n s V(p = 2(fNn s 
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As the result we find an equation where the contributions proportional to ip cancel each 
other. 



d t ip{x,t) = D(x)V 2 ^(x,t) 
with D(x) = {z- l)n s (x) (16) 



As before let us study a radial symmetric solution. Taking into account eq. (|bl| ) and inspect- 
ing the last equation in comparison with the diffusion equation we make the ansatz 



r a+2 



oc t~"< exp(-b-) (17) 

which allows to determine the unknown quantities 7 and b: 

a + d 1 

7 = and b= — — — (18) 

' a + 2 c(z-l){a + 2) 2 v ; 

Like before the result leads to the conventional behaviour for d = 2. In the limit of an 
infinite range of the hopping process corresponding to z — > 00 we find 7 = 2. 
Summarizing the results the radial symmetric solution of eq.flTUD can be written for large 
distances r in scaling form 

n (r,t)~±+r'/(^) (19) 

where the exponents and the scaling function are given by 

p= Z -^\ f(v) = Av {2a)na+2)iz - 1) exp(-to) (20) 
a + 2 

To derive the last equation we have also used the radial symmetric solution of the function 
0(r) introduced after eq. (|T3 ) . 

Inspecting eq. C|XQ|) one remarks that the dynamical part of the solution yields a maximum. 
The position of this bump r& follows a simple law 

r b {t) oc (21) 

Obviously, the bump is the more pronounced the lower the dimension. The result is depicted 
in Fig.l. 

The exponent of diffusion d w is simply defined by eq. (|21|) 

d w = a + 2 (22) 

As mentioned the exponent agrees with that of conventional diffusion only in the two dimen- 
sional case. For instance for d = 3 it results d w = We anticipate that d w is considerably 
larger than 2 since many of the sites are already occupied due to the restriction introduced 
by eq.dH). As the consequence a lot of the neighbouring sites of a given pair of sites are 
unavailable to be reached by particles. Instead of that they are obliged to return to the 
starting point. The restrictions act like a cage for the particles. A similar conclusion can 
be drawn in case of a fractal lattice which leads also to an exponent d w larger than two. 
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The position of the bump r& moves much less rapidly with time t than for an unrestricted 
hopping. 



In the next step we would like to discuss the influence of a heat bath with a temperature 
T. To this aim the Fock space formulation is extended by replacing the evolution operator 
Lby§ 

L = Jij [(1 - did]) exp(-H/2T)4d j exp(F/2T)] (23) 

<ij> 

where the Hamiltonian will be chosen to be 

H = hJ2 n r + K J2 n r n s (24) 

r <r.s> 

Here h means an external field (or chemical potential of the lattice gas) and K stands for 
the static coupling between nearest neighbours situated at the cells r and s. For K > 
there is an additional repulsion whereas for K < the particles are attractive. Remark that 
eq. (p3|) is in accordance with detailed balance. 

Using the algebraic properties of Pauli-operators the evolution operator can be rewritten 

L — (1 — did])d\djAiBj 

<ij> 

with Ai = exp I -iT/T^n, I and Bi = A' 1 (25) 
V »W / 

The exact evolution equation for the density reads now 

d t (n r ) = exp(K/T) $3<J P j(E/rV 2 VS - V r V 2 Uj)) 

with U r = (1 — n r )A r and VJ. = n r B r (26) 

As far as the hopping parameter does not depend on the concentration itself the last equa- 
tion can be represented as a continuity equation. However, in the case where the Jy are 
determined by the concentration cf. eq.(§) the conservation law is not valid as already re- 
marked in the non-thermalized approach. After performing the continuum limit we obtain 
in the one dimensional case 

d t n = nd xx n - 2(d x n) 2 + A{K/T){nd xx n - (d x n) 2 - n 2 d xx n) (27) 

Applying the same procedure as before the last equation can be transformed in 

d t ip = D(x)d xx ijj 

with D(x) =n s (x)(l + A(K/T)n s (x)(l~n s (x)) (28) 

The only difference to eq. (TO) is that D(x) is replaced by D(x). If K > 0, i.e. repulsive in- 
teraction, the modified diffusion function D(x) is always positive and the general behaviour 
is only slightly changed. However if K < (attraction) the quantity D(x) can be negative 
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depending on the ratio ^ and on the density n s . A negative diffusion coefficient is obviously 
related to an instability. 

In the present paper we have demonstrated that constraints during the diffusion process lead 
to a complete different behaviour. Due to the restrictions imposed on the diffusion process 
it appears a non-linear diffusion equation which can be analyzed partially in an analytical 
manner. There appears a maximum in the concentration which follows anomalous diffusion. 
Such a behaviour should be relevant for glassy materials characterized by rapid decrease of 
the mobility for decreasing temperatures. Further applications can be discussed when other 
degrees of freedom as impurities or vacancies are taken into account. 
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FIGURES 
FIG. 1. 3 dimensional plot n(r, t) for d = 1 
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